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Abstract 



We discuss a few selected topics in finite temperature field theory. 
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1 Introduction 



Studies of physical systems at finite temperature have led, in the past, to many interesting proper- 
ties such as phase transitions, blackbody radiation etc. However, the study of complicated quantum 
mechanical systems at finite temperature has had a systematic development only in the past few 
decades. There are now well developed and well understood formalisms to describe finite tem- 
perature field theories, as they are called. In fact, as we know now, there are three distinct, but 
equivalent formalisms [1-3] to describe such theories and each has its advantages and disadvantages. 
But, the important point to note is that we now have a systematic method of calculating thermal 
averages perturbatively in any quantum field theory. 

This, of course, has led to a renewed interest in the study of finite temperature field theories 
for a variety of reasons. We can now study questions such as phase transitions involving symmetry 
restoration in theories with spontaneously broken symmetry [4]. We can study the evolution of 
the universe at early times which clearly is a system at high temperature. More recently, even 
questions such as the chiral symmetry breaking phase transition or the confinement-deconfinement 
phase transition in QCD [5-6] have drawn a lot of attention in view of the planned experiments 
involving heavy ion collisions. This would help us understand properties of the quark-gluon plasma 
better. 

The goal of this article is to share, with the readers, some of the developments in finite temper- 
ature field theories in the recent past and the plan of the article is as follows. In the next section, 
we will describe some basic ideas behind describing a quantum mechanical theory in terms of path 
integrals [7]. This is the approach which generalizes readily to the study of finite temperature field 
theory. In section 3, we will discuss one of the formalisms, in fact, the oldest one, of describing 
finite temperature field theory. This goes under the name of the imaginary time formalism or the 
Matsubara formalism [1, 5, 8-10]. In this description, the dynamical time is traded in for the tem- 
perature. In contrast, the real time formalisms of finite temperature field theory contain both time 
and temperature. In section 4, we discuss one of the real time formalisms known as thermo field 
dynamics [3, 10-11]. This is an ideal description to understand operator related issues involving 
finite temperature field theories although it has a path integral representation which is quite nice 
for calculations as well. The other real time formalism, which is much older and is known as the 
closed time path formalism [2, 10, 12], is described in section 5. This formalism is very nice because 
it describes both equilibrium and non-equilibrium phenomena, at finite temperature, with equal 
ease. Temperature leads to many subtle features in field theories. In section 6, we discuss one such 
subtlety, namely, how one needs a generalization of the Feynman combination formula to perform 
calculations at finite temperature [13]. In section 7, the issue of large gauge invariance is discussed 
within the context of a simple quantum mechanical model [14-15]. In section 8, we discuss in some 
detail how temperature can lead to breaking of some symmetries like supersymmetry [16] (Tem- 
perature normally has the effect of restoring symmetries). Finally, we present a brief conclusion 
in section 9. The subject of finite temperature field theories is quite technical and to keep the 
contents simple, we have chosen, wherever possible, simple, quantum mechanical models to bring 
out the relevant ideas. Finally, we would like to note that there are many works in the literature 
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and the references, at the end, are only representative and are not meant to be exhaustive in any 
way. 



2 Path Integrals at Zero Temperature 

In studying a quantum mechanical system or a system described by a quantum field theory, we 
are basically interested in determining the time evolution operator. In the standard framework of 
quantum mechanics, one solves the Schrodinger equation to determine the energy eigenvalues and 
eigenstates simply because the time evolution operator is related to the Hamiltonian. There is an 
alternate method for evaluating the matrix elements of the time evolution operator which is useful 
in studying extremely complicated physical systems. This goes under the name of path integral 
formalism [7, 17-18]. 

In stead of trying to develop the ideas of the path integral formalism here, let us simply note 
that, for a bosonic system described by a time independent quantum mechanical Hamiltonian, the 
transition amplitude can be represented as (The subscript H denotes the Heisenberg picture.) 

H{xf,tf\x,,U)H = (x/|e-^^(*/-*>)|xi) = J Pxe^^M (1) 

There are several comments in order. First, the transition amplitude is nothing other than the 
matrix element of the time evolution operator in the coordinate basis. Second, the integral on the 
right hand side is known as a path integral. It is an integral over all possible paths connecting the 
initial coordinate Xj and the final coordinate xj which are held fixed. The simplest way to evaluate 
such an integral is to divide the time interval of the path between Xi and x j into intervals of equal 
length. Integrating over all possible values of the coordinates of the intermediate points (which 
are ordinary integrals) and taking N ^ oo such that the time interval is held fixed is equivalent 
to integrating over all possible paths. Finally, the action S[x] in the exponent of the integrand is 
nothing other than the classical action for the bosonic system under study. This is true for most 
conventional physical systems where the Hamiltonian depends quadratically on the momentum. If 
this is not the case (and there are some cases where it is not), the right hand side of (Q) needs to 
be modified. However, for most systems that we will discuss, we do not have to worry about this 
fine point. 

The advantage of the path integral is that while the left hand side involves quantum mechanical 
operators, the right hand side is described only in terms of classical variables and, therefore, the 
manipulations become quite trivial. Furthermore, the transition amplitude defined in eq. (|I|) can 
be generalized easily to incorporate sources and this allows us to derive various Greens functions 
of the theory in a very simple and straightforward manner. As an example, let us simply note here 
that for a harmonic oscillator, the action is quadratic in the dynamical variables, namely, 
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and, in this case, the path integral can be exactly evaluated and has the form [7] 



muo 



lirih sinwT 



(2) 



Here, we have defined T = tf—ti. S[xci] represents the action associated with the classical trajectory 
(satisfying the Euler-Lagrange equation) and has the form 



Six. 



muj 



cll 



2 sin ujT 



{xf + x^f) cos ujT — 2xiX f 



(3) 



The path integrals can also be extended to quantum mechanical systems describing fermionic 
particles. However, one immediately recognizes that there are no classical variables which are 
fermionic. Therefore, in order to have a path integral description of such systems in terms of 
classical variables, we must supplement our usual notions of classical variables with anti-commuting 
Grassmann variables [19]. With this, for example, we can write a classical action for the fermionic 
oscillator as ^ 

S[iP, V^] = f ^ dt {iijip - ujipip) (4) 
Jti 

Here ip and ^ are anti-commuting Grassmann variables and in the quantum theory, as operators, 
can be identified with the fermionic annihilation and creation operators respectively. The action in 
eq. (^) is also quadratic in the variables much like the bosonic oscillator and the path integral for 
the fermionic oscillator can also be exactly evaluated giving [7] 



e 2 



In a quantum field theory, we are often interested in evaluating time ordered correlation func- 
tions in the vacuum because the S-matrix elements can be obtained from such Greens functions. 
These can be derived in a natural manner from what is known as the vacuum to vacuum transition 
functional which can be obtained from the transition amplitude in eq. (|l]) in a simple manner and 
also has a path integral representation of the form 

^lim (0|e-^^^|0) = j Pxet^M (6) 

where 

/oo 
dtL{x,x) (7) 
-oo 

Furthermore, the path integral in eq. (^) has no end-point restriction unlike in eq. (|l]). This 
vacuum to vacuum transition amplitude is also commonly denoted by (0|0) with the limiting process 
understood. We note here that an analogous formula also holds for fermionic systems. 
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The vacuum to vacuum amplitude in the presence of a source has the form 

Z[J] = (0|0)j = J Dxes^l^'-^l (8) 

where 

/oo 
dt J{t)x{t) (9) 
-oo 

Here J{t) is a classical source and it can be easily checked that, in the limit of vanishing source, 
the functional derivatives of Z[J] give rise to time ordered Greens functions in the vacuum. 

With this very brief review of the path integral description for zero temperature quantum me- 
chanical theories, we are now ready to describe the different formalisms available to study quantum 
mechanical systems at finite temperature. 



3 Imaginary Time Formalism 



The properties of a quantum mechanical system, at finite temperature, can also be given a path 
integral description. There are various, but equivalent ways of doing this. Of the different for- 
malisms available to study a quantum mechanical system at finite temperature, the imaginary time 
formalism is the oldest [1]. To appreciate this, let us recall some of the features of a statistical 
ensemble. A statistical ensemble in equilibrium at a finite temperature ^ (in units of Boltzmann 
constant) is described in terms of a partition function 

Z{I3) = TV = TV e-^'^ (10) 

Here p{j3) is known as the density matrix (operator) and Ti can be thought of as the generalized 
Hamiltonian of the system. If 

rL = H 

where H is the Hamiltonian of the system, we say that the ensemble is a canonical ensemble where 
the particle number is fixed and the system is allowed to exchange only energy with a heat bath. 
On the other hand, if 

where N is the number operator, then, the ensemble is known as a grand canonical ensemble where 
the system can exchange not only energy with a heat bath, but can also exchange particles with a 
reservoir. The constant fi is known as the chemical potential. In a statistical ensemble, of course, 
the important observables are the ensemble averages and, for any observable O, they are defined 
as 

(C?)/3 = ^^P(/3)C? (11) 

Let us also note here that since the partition function involves a trace, it leads to an interesting 
identity following from the cyclicity of the trace, namely, (we will assume from now on, unless 
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otherwise specified, that h = 1) 



{0,{t)02{t'))p = J^Tre-^«Oi(t)O2(0 



m 

= -^^Tr e-^'^02{t')Oi{t + il3) 

= {02{t')0i{t + iP))p (12) 

Such a relation is known as the KMS (Kubo-Martin-Schwinger) [20] relation which generalizes to 
all statistical ensemble averages and plays a crucial role in the study of finite temperature field 
theories. 

It was observed quite early by Bloch [21] that the operator e^^^ in the definition of the partition 
function is like the time evolution operator in the imaginary time axis. This is really at the heart 
of the imaginary time formalism. In fact, let us note that the canonical partition function can be 
written as (with the trace taken in the coordinate basis) 

Z(/3) = / dx{x\e-'^^\x) (13) 



It is clear now that if we identify T = —i(3 in eq. (2), then, we can give the partition function a 
path integral representation as {h = 1) 

Z{P) = [ Pxe-^^t^l (14) 



where 5'£;[x] is the Euclidean (imaginary time) action for the system defined over a finite time 
interval as 

Se[x] = / dtLE{x,x) (15) 

^0 



Furthermore, it is clear from eq. (13) that the variable x must satisfy the periodic boundary 
condition 

x{P) = x(0) (16) 

for eq. (|l4|) to represent a trace (namely, the initial and the final states must be the same) and 
that the end point is being integrated over in the path integral in eq. (|l^ unlike in eq. (||). (It is 
important to note that the original work of Matsubara is an operator description of the imaginary 
time, but we will not discuss it in the present article.) 



In fact, as an example, let us evaluate the canonical partition function for the bosonic oscillator 
using this formalism [7]. The transition amplitude is already given for zero temperature in eq. (^). 
Now making the identifications 

T = -i/3, Xi = Xf = x (17) 
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we obtain from eqs. (0) and (13) 



g— (mu) tanh ^ ) 



1 / 

mw \ 2 / vr 



27r sinh /3u; / y muj tanh ^ / 

e 2 

(18) 



e/3(^ - 1 

This is, indeed, the partition function for the bosonic osciUator as can be directly verified 



The partition function, for a fermionic system, can also be similarly given a path integral 
representation. However, the anti-commuting nature of the fermion variables introduces one crucial 
difference, namely, for a fermion theory, we have 

Z{(3) = J VtpVi^e-^^^^'^^ (19) 

with anti-periodic boundary conditions [10] 

^PiP) = -m, = -m (20) 

The Euclidean (imaginary time) action is again defined over a finite time interval as in eq. (^). 
In fact, let us calculate the canonical partition function for a fermionic oscillator, as an example, 
from the result in eq. (|5|) as well as the identifications in ( [20| ) [7]. Using 

we obtain (remember T = —if3) 

= (l + e"^'^) =2cosh^ (21) 

In evaluating this, we have made use of the Berezin rules of integration [19] for Grassmann variables 
and we note that eq. (^l]), indeed, gives the correct partition function for a fermionic oscillator as 
can be directly calculated. 



Although our discussion so far has been within the context of simple quantum mechanical 
systems, everything we have said can be carried over to a quantum field theory. The partition 
function for a quantum field theory can again be written as a path integral involving a Euclidean 
action as 

Z{P) = J V^Vi)V(t)e-^^^^^^^^^ (22) 

where the Euclidean action is defined over a finite time interval and the fields satisfy the periodicity 
(anti-periodicity) conditions 

4>{I3, x) = </.(0, x), V(/?, S) = -V'(0, x) (23) 
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and so on. The discussion is slightly more involved for gauge theories and to keep things simple, 
we will not discuss gauge theories. 



This formulation of a field theory at finite temperature is known as the imaginary time formalism 
or the Matsubara formalism [1] and is the oldest formalism. There are several distinguishing features 
of this formalism. For example, since the time interval is finite, Fourier transformation of the time 
variable would involve discrete energies. In other words, the Fourier transform of the propagator, 
say for example, at finite temperature in the imaginary time formalism, would take the general 
form 

^/3(^,^) = 4Ee"*"""^/5K'^) (24) 

P n 

where u;„ = with n = 0, ±1, ±2, • • •. However, from the definition of the time ordered product 
r.(</.(r)</.t(T')) = Qir - t')4>{t)4>\t') ± dir' - t)4>\t')<^{t) (25) 



where we have allowed for both bosonic and fermionic fields and the KMS condition in eq. (|T2|), it 
follows that, for r < 0, 

g/3(T,f) =±e?/3(T + /3,f)) (26) 

It is important to recognize that the periodicity (anti-periodicity) of the propagator arises from the 
definition of the time ordered product for the bosonic (fermionic) fields and the KMS condition and 
is not directly connected with the periodicity (anti-periodicity) of the corresponding field variables 
which we have discussed earlier. This periodicity (anti-periodicity) of the propagator, on the other 
hand, leads to the restriction that eq. (|23) holds with 



for bosons 

where n = 0, ±1, • • •. These are conventionally known as the Matsubara frequencies [22]. 



n \ (2ra+i)7r fermions ^ '^^ 



Given this, one can now calculate the propagators for bosonic and fermionic field theories in 
the Matsubara formalism and they take the forms (in the momentum space) 

QfiiyJnS) = — i ; = ,2n7rN9 ^9 7 (^8) 

c, A - 7V + 7-fc + r. _ /(^)+7-fe + m 



Perturbative calculations can now be developed quite analogously to the zero temperature field 
theory. For example, given a field theory, we can read out the vertices from the Euclidean form of 
the action and use the propagators of eq. (^, [29| ) to carry out a diagrammatic calculation which 
would lead to the ensemble average for a given observable. It is clear that, because the time interval 
is finite in this formalism, the coordinate space calculation of any diagram is cumbersome. However, 
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much like at zero temperature, the momentum space calculation is much simpler. However, one 
should keep the difference in mind, namely, that, at finite temperature, the external and the internal 
energies are discrete as in eq. (p7|). Consequently, the integration over internal energies (of zero 
temperature) is replaced by a sum over the internal energies. More specifically, we must use 




As an example, let us consider the self-interacting scalar theory described by 

We note that the only one loop correction in this theory is the mass correction. Rotating to 
Euclidean space and using the propagator for a scalar theory as given in eq. (|^) as well as (|^), 
we obtain the one loop mass correction to be 

2/3^7 (2^)3 (2|ZL)2 + ^2+^2 

= V^VV/^'^ ^ (32) 

2/3 (2vr)3 n2 + (^)2 ^ ' 

Here, we have introduced the notation, 

u;fc = (P+m2)5 (33) 
The sum, in eq. (p2[), can be easily evaluated using the method of residues leading to 

y 7T = — coth Try for y > (34) 



n=—oo 



Using this, the one loop mass correction can be determined to be [23] 



Am = — - — -TT — coth 



4 J (27r)3 oJk 

+ 



4 J (27r)3 wfc 2 J (27r)3 uJk 6^'^'= - 1 

2 

/3 



Ami + ^"^fl (35) 



There are several things to note from this calculation. First, the mass correction separates 
into two parts - one independent of temperature and the other genuinely a finite temperature 
correction. The temperature independent part (zero temperature part) is divergent as is expected 
at zero temperature and the divergence has to be handled by the usual process of renormalization. 
However, the finite temperature part is completely free from ultraviolet divergence. This is a general 
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feature of finite temperature field theories that temperature does not introduce any new ultraviolet 
divergence. We will return to this question later within the context of real time formalisms for finite 
temperature field theories. Let us also note that (see (|35D) the finite temperature integrals are, in 
general hard to evaluate and cannot be evaluated in a closed form. However, we can always make 
a high temperature expansion (small (3) which would give the temperature dependent correction to 
the mass as 

" 2^55 = ^ (36) 

This shows that temperature induces a mass correction which is positive. Intuitively, it is clear 
that this is the behavior we would expect from a particle moving in a medium and, furthermore, 
the positivity of this correction is crucial in the study of symmetry restoration in field theories with 
spontaneous symmetry breaking. 



This gives a flavor of calculations at finite temperature, particularly, in the imaginary time 
(Matsubara) formalism. It is worth noting here that, by construction, the imaginary time formalism 
would describe physical systems in equilibrium quite well. Since we have traded the time variable 
for temperature, it is well suited to calculate static, equilibrium quantities. Slow temperature 
dependence can, however, be brought in by analytically rotating the final result to Minkowski time 
[24]. This rotation is, on the other hand, nontrivial since we only have information about quantities 
at discrete energy values in the Euclidean space. The imaginary time formalism is not suitable to 
discuss non-equilibrium phenomena. 



4 Thermo Field Dynamics 



As we have seen, in the imaginary time formalism, the time variable is traded for the temperature. 
However, in studying various processes, it is desirable to have the time coordinate in addition to 
the temperature. Formalisms where this can be achieved are known as the real time formalisms 
and there are two distinct, but equivalent such formalisms. In this section, we will discuss the 
formalism of thermo field dynamics [3, 11, 25] returning to the alternate formalism in the next 
section. 

Let us recall from ([l^) that the ensemble average of any observable is given by 

= 2^1: ("I0|n> (37) 

Here, we have assumed that the eigenvalues of Ti. are discrete, for simplicity, and that 

H\n) = En\n) 
{m\n) = 6mn 
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^|n)(n| = I 

n 



(38) 



At zero temperature, we know that the Feynman diagrams correspond to vacuum expectation 
values of time ordered products. Thus, intuitively, it is clear that if we can express the ensemble 
averages as expectation values in some vacuum (say, a thermal vacuum), then, we can take over 
all the diagrammatic machinery of the zero temperature field theory. The question, therefore, is 
whether we can define a vacuum, say |0,/3), such that we can write any ensemble average as 

{0)p = (0,/3|O|0,/3) = -i- ^e-^^" {n\0\n) (39) 



Let us suppose that we can define such a thermal vacuum state as a linear superposition of the 
states in our physical Hilbert space, namely, 

|0,/3) =5]|n)(n|0,/3) =^/„(/3)|n) (40) 

n n 

This would lead to 

(0,/3|O|0,/3) /:(/?)/„(/?) (n|0|m) (41) 

n,m 

Consequently, this would coincide with eq. (|39|) only if 

fmfm{f3) = ^^e-P^-5mn (42) 

Since /„'s are ordinary numbers and eq. (^) is more like an orthonormality condition, it is clear 
that we cannot satisfy this condition (and, therefore, define a thermal vacuum with the right 
properties) if we restrict ourselves to the original Hilbert space. 



On the other hand, it is also clear from this analysis that if /„'s, somehow, behave like a state 
vector, then, the condition in eq. (^) can be easily satisfied. In fact, let us introduce a fictitious 
system identical to our original system and denote it by a tilde system. The states in the combined 
Hilbert space of this doubled system would have the form 

|n, fh) = \n) \ fn) 

Let us assume that the thermal vacuum can be written as a linear superposition of states in this 
doubled Hilbert space of the form 

|0, /5> = E = E ® 1^) (43) 

n n 

This would lead to 

(0,/3|O|0,/3) = Y.fniP)fmm{n,n\0\m,m) 

n,m 

= J2fni(3)fM{n\0\m)Sn,m 

= Y.fnmn{f3){n\0\n) (44) 

n 
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In deriving this result, we have used the fact that an operator of the original system does not 
act on states of the tilde system and vice versa. The result in eq. (44) is quite interesting because 
it says that if we choose 



f,-l3E„ p-l3E„/2 

fmfn{f3) = ^^ or, /„(/?) = /:(/?) = (45) 
then, eq. ( |4^ ) would, indeed, coincide with the ensemble average in eq. (|39[). 

This analysis shows that it is possible to introduce a thermal vacuum such that the ensemble 
average of any operator can be written as the expectation value of the operator in the thermal 
vacuum. The price one has to pay is that the Hilbert space needs to be doubled. The advantage, 
on the other hand, lies in the fact that the description would now involve both time and temper- 
ature (since we have not traded time for temperature) and all the diagrammatic methods of zero 
temperature field theory can now be taken over directly. 



Fermionic Oscillator 



To get a flavor for things in this formalism, let us analyze in some detail the simple quantum 
mechanical system of the fermionic oscillator. The Hamiltonian for the system is given by {fi = 1) 



H = uja)a 



(46) 



Here, the fermionic creation and annihilation operators satisfy the canonical anti-commutation 
relations 



a, a) 
[a, a] 



1 



a\ a) 







(47) 



In this case, the spectrum of the Hamiltonian is quite simple and the Hilbert space is two dimen- 
sional with the basis states given by |0) and |1) = a^^jO). 



According to the general philosophy of thermo field dynamics, we are supposed to introduce a 
fictitious tilde system which is identical to our original system. Thus, we define 



H = ujo^a 



(48) 



with the anti-commutation relations 



a, a' 



a\ a) 



(49) 



12 



Furthermore, we assume that the creation and the annihilation operators for the tilde and the 
non-tilde systems anti-commute. 



The Hilbert space for the combined space is now four dimensional and, following our earlier 
discussion, we choose the thermal vacuum to be 

\o,(3) = Wm ^ |6) + ® |i) (50) 

The normalization of the thermal vacuum gives 

(0,/3|0,/3) = |/o(/3)P + |/i(/3)|2 = l (51) 
while the expectation value of the number operator gives 

(0,/3|iV|0,/3) = (0,/3|ata|0,/3) = \Mf3)\^ = (52) 

From these, we can obtain 

SO that we can write 

To further understand the properties of this system, let us note that we can define a Hermitian 
operator in this doubled space 

G{e) = -ie{j3) {aa - a) a}) (55) 
This would, in turn, lead to a formally unitary operator 

U{P) = e-^^(^) (56) 
which would connect the thermal vacuum to the vacuum of the doubled space, namely, 

C/(/3)|0,0) = cose(/5)|0,0) + sin^(/3)|l,i) = |0,/3) (57) 

provided 

cos0(/?) = = / sin0(/3) = /i(/?) = ^ " ' (58) 

The unitary operator would also induce a transformation on the operators of the form 

0{P) = U{P)OU\l3) (59) 

In particular, this would give 

a{P) = cos 9 {(3) a -sin 9 {P)a^ 

a{P) = cos9{(3)a + sm9{P)a^ (60) 
13 



as well as their Hermitian conjugates. These operators would satisfy the same anti-commutation 
relations as the original ones and we can think of them as the thermal creation and annihilation 
operators. Consequently, we can build up the thermal Hilbert space starting from \0,(3) and the 
thermal creation operators. 

In particular, it is trivial to check, using (^) that the thermal vacuum satisfies 
a(/3)|0,/3) = {cose{/3)a-sm9{p)d^)\0,p) =0 

a(/3)|0,/3) = (cos0(/?)5 + sin0(/?)at)|O,/3) = (61) 

This is quite interesting for it says that annihilating a particle in the thermal vacuum is equivalent 
to creating a tilde particle and vice versa. Consequently, we can intuitively think of the tilde 
particles as kind of hole states of the particles or particle states of the heat bath. This gives a nice 
intuitive meaning to the doubling of the degrees of freedom in thermo field dynamics. Namely, an 
isolated system in thermal equilibrium really consists of two components - the original system and 
the heat bath. 



We also note here that although the operator connecting the thermal vacuum to the vacuum 
in the doubled space is formally unitary, it is more like a Bogoliubov transformation. In more 
complicated models with an infinite number of degrees of freedom (namely, in field theories) such 
an operator takes us to a unitarily inequivalent Hilbert space. Let us also note here, for future use, 
the simple formula following from eq. ( pO| ) that 

= U{P) ( \ (62) 




where 

n(R\ - I 

sine(/3) cose{(5) 



Finally, let us conclude the discussion of this example by noting that the states in the thermal 
Hilbert space are eigenstates of neither H nor H. Rather they are the eigenstates of the operator 

H = H- H (64) 

Furthermore, this combination of the Hamiltonians is also invariant under the unitary transforma- 
tion of (|57|). This is, indeed, the Hamiltonian that governs the dynamics of the combined system. 



Bosonic Oscillator 



The analysis for the case of the bosonic oscillator is quite analogous to the discussion of the fermionic 
oscillator. Therefore, without going into too much detail, let us summarize the results. First, the 
Hamiltonian for the system is given by 

H = ua^a (65) 
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much like the fermionic oscillator. However, the creation and annihilation operators satisfy canon- 
ical commutation relations of the form 



a, a) 
[a, a 



a\ a) 







(66) 



The Hilbert space for the bosonic oscillator is infinite dimensional with the energy eigenstates given 
by 

H\n)=nu)\n), re = 0,1,2,- •• (67) 



According to the general discussions of thermo field dynamics, we introduce an identical, but 
fictitious tilde system with the Hamiltonian 



H = ujo^a 



(68) 



The tilde creation and annihilation operators are expected to satisfy commutation relations analo- 
gous to (|66|). Furthermore, the tilde operators are supposed to commute with the original operators 
of the theory. 



Following the discussion of the earlier section, we can determine the thermal vacuum state in 
this case to be 

oo 

|0, /3) = (1 - e-^'")l/2 ^ ^gg) 
n=0 

As in the fermionic oscillator, we can introduce the Hermitian operator 

G{9) = -ie{l3) (aa - a^a)) (70) 

and the unitary operator 

[/(/?) = e-*^W (71) 

Then, it is straightforward to check and see that the unitary operator connects the thermal vacuum 
to the vacuum of the doubled space provided 



The unitary operator induces a transformation of the operators of the form 

0(/3) = C/(/3)OC/t(/?) (73) 

leading to 

a{l3) = cosh9{(3)a-sinhe{l5)d) 
a{l3) = coshe(/3)o-sinh0(/3)a"f 



15 



and similarly for the Hermitian conjugates. As we have seen in the last section, these can be 
thought of as the creation and annihilation operators for the thermal Hilbert space. In particular, 
the thermal vacuum is easily seen to satisfy 

o(/3)|0,/3) = {coshe{/3)a-sinhe{/3)a^)\0,(3) =0 

a{p)\0,P) = (cosh0(/3)a-sinh0(/3)at)|O,/3) =0 (74) 

This, again, reinforces the intuitive picture of doubling in thermo field dynamics. Let us also note 
here, for future use, the simple formula following from ( [7^ ) 

( Z ) - ( r. ) p=) 

where 

ff(n^_( cosh6'(/?) -smh9{P)\ , 

^^^'^ \ -sinh^(/?) cosh0(/3) ) 



Field Theory 



The extension of these results to a field theory is quite straightforward once we keep in mind that, 
at the free level, a quantum field theory is simply an infinite collection of oscillators with frequencies 
dependent on the momentum of the mode. Consequently, the thermal vacuum, in this case, would 
be connected to the vacuum of the doubled space as 

|0,/3) = C/(/3)|0,0) = e-'^(^) 1 0,0) (77) 

where 

Gi(^) = E hiP) [hH - (78) 

k 

with, say, for bosons, 

coshOrW) = , sinherr(/3) = , (79) 

Here, for a relativistic theory, we have 



Let us next note that, at zero temperature, the original fields are decoupled from the tilde fields. 
Thus, if we were to define a doublet of fields (real scalar field) as in eq. ([f^) 

(80) 
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then, at zero temperature the propagator is defined to be (This is not to be confused with the 
generator of the BogoUubov transformations in eqs. (|55|) , (|70|) and ([78|)) 



iG{x-y) = (0,0|r(cl.(x)cl>(y))|0,0) 
which has the momentum space representation 

\ k'-^—m'-^—it / 

Given this, the finite temperature propagator can be determined to be 

iGp{x-y) = (0,/3|r($(:E)$(y))|0,/3) 

= {0,0\uHP)TmxMy))U{(3)\0,6) (82) 



Using now the generahzation of eqs. ( |73| , [75| , [f6D , the momentum representation for the propagator 
can be determined to be 

Gp{k) = u{-e^)G{k)u^{-e^) 

= (^^^ \ )-2zvrnM|fe°|)^(fe^-m^)( ,|,0|i ) (83) 



There are several things to note from the structure of the propagator which are quite general 
for a real time formalism. First, the propagator is a 2 x 2 matrix, a consequence of the doubling 
of the degrees of freedom. Second, the propagator is a sum of two parts - one representing the 
zero temperature part and the other representing the true temperature dependent corrections. 
The propagator is still the Greens function for the free operator of the theory, but corresponding 
to different boundary conditions (remember the KMS condition in eq. ([T^)). While the zero 
temperature part of the propagator corresponds, as usual, to the exchange of a virtual particle, the 
temperature dependent part represents an on-shell contribution (because of the delta function). In 
fact, the intuitive meaning of the temperature dependent correction is quite clear. In a hot medium, 
there is a distribution of real particles and the temperature dependent part merely represents the 
possibility that a particle, in addition to having virtual exchanges, can also emit or absorb a real 
particle of the medium. 

Since the temperature dependent part of the propagator is on-shell, it is clear that there can 
be no new ultraviolet divergence generated at finite temperature. All the counter terms needed 
to renormalize the theory at zero temperature would be sufficient for studies at finite temperature 
as well. (Of course, the infrared behavior is another story. Infrared divergence, in a field theory, 
becomes much more severe at finite temperature, a topic that I will not get into.) There is an alter- 
nate way to visualize this. At finite temperature the distribution of the real particles is Boltzmann 
suppressed as we go up in energy and, consequently, thermal corrections corresponding to infinite 
energy cannot arise. 

Once, we have the propagator, we can venture to do a diagrammatic calculation in this for- 
malism. The only things missing are the interaction vertices of the theory. There is a well defined 
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procedure [26] (called the tilde conjugation rule) to construct the complete Lagrangian from which 
to construct the vertices. Very simply, it corresponds to what we have noted earlier, namely, the 
dynamical Hamiltonian (and, similarly, the Lagrangian) is as given in eq. (|64[). It is simply the 
difference between the original and the tilde Hamiltonians. Thus, we see that the complete theory 
would contain two kinds of vertices - one for the original fields while the second for the tilde fields. 
The vertices for the tilde fields will have a relative negative sign corresponding to the original 
vertices. Given the vertices and the propagator, it is now straightforward to carry out any dia- 
grammatic calculation to any order. Let me emphasize here that although, at the tree level, there 
is no vertex containing both the original and the tilde fields, such vertices would be generated at 
higher loops because of the nontrivial matrix structure of the propagator. 

Thermo Field dynamics is a real time formalism. But, more than that, it is really an operator 
formalism and hence very well suited to study various operator questions such as the structure 
of the thermal vacuum, the theorems on symmetry breaking etc. It can also be given a path 
integral representation and corresponds to choosing a specific time contour in the complex t plane 
[27, 10] (remember that the imaginary time formalism also corresponds to choosing a specific time 
contour, namely, along the imaginary time axis) and I will come back to this question in the next 
section. However, once again from the philosophy of thermo field dynamics, it is clear that, it 
is a natural formalism to describe equilibrium phenomena where quantities depend on both time 
and temperature. While there are several attempts to generalize this to include non-equilibrium 
phenomena, there does not yet exist a complete description. 



5 Closed Time Path Formalism 



The closed time path formalism is also a real time formalism which was formulated much earlier than 
thermo field dynamics within the context of non-equilibrium phenomena [2]. The two formalisms 
are, in some sense, complementary to each other although the closed time path formalism can 
describe both equilibrium and non-equilibrium phenomena with equal ease. 

The basic idea behind the closed time path formalism [10] is the fact that when a quantum 
mechanical system is in a mixed state, as is the case in the presence of a heat bath, the system can 
be naturally described in terms of a density matrix defined, in the Schrodinger picture, as 

pit)=J2Pn\Mt)){Mt)\ (84) 

n 

Here, pn represents the probability for finding the quantum mechanical system in the state |^„(t)) 
and, for simplicity, we have assumed the quantum mechanical states to form a discrete set. It is 
Pn which contains information regarding the surrounding which is hard to determine, but, being a 
probability, it satisfies 

n 
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Given the density matrix, the ensemble average of any operator can be calculated in the 
Schrodinger picture as 

{0){t) = J2pn {Mt)\0\Mt)) = Trp(t) O (85) 

n 

The ensemble average, in this case, naturally develops a time dependence from the time dependence 
of the density matrix. In this formalism, we can naturally define an entropy as 

S = - '^Pn Inpn 

n 

which is by definition positive semi-definite and measures the order (or lack of it) in an ensemble. 



The state vectors satisfy the Schrodinger equation {h = 1) 

From this, we can determine the time evolution of the density matrix which turns out to be the 
Liouville equation 

^^ = [H,p{t)] (86) 

In deriving this, we have assumed that the probabilities do not change with time (appreciably) 
implying that entropy remains constant during such an evolution. The reason for this assumption 
is our lack of knowledge about the time evolution of the surrounding such as the heat bath. On 
the other hand, adiabatic evolutions do arise frequently in physical systems and, consequently, we 
would continue with this assumption. 



Let us note that eq. ( |86[ ) has a simple solution of the form 

p{t) = U{t, 0)p(0)[/t {t, 0) = U{t, 0)p{0)U{0, t) (87) 
where the time evolution operator has the general form 

C/(f,t')=T(^e-^/>*"^(*")) (88) 

Furthermore, it satisfies the semi-group properties 

U{ti,t2)U{t2M) = 1 

U{ti,t2)U{t2M) = U{ti,h) ioTti>t2>h (89) 

In particular, let us note that if the Hamiltonian is time independent, eq. (^^ takes the simple 
form 

p{t) = e-'"'p{0)e'^' 

and, furthermore, if the Hamiltonian commutes with p{0), the density matrix would be time in- 
dependent, describing a system in equilibrium. This would be true, for example, if the states in 
eq. ( |8^ are stationary states. This is also true if the probabilities have a Boltzmann distribution 
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in which case, we refer to the system as being in thermal equihbrium. However, we will not re- 
strict to any such special case allowing for the formalism to accommodate both equilibrium and 
non-equilibrium phenomena. 

Keeping in mind the fact that we are ultimately interested in a thermal ensemble, let us choose 

for some Hi. Since the density matrix is a positive Hermitian matrix with unit trace, mathemat- 
ically, this is allowed. But, more important is the physical reason behind such a choice. Namely, 
we can think of the dynamical Hamiltonian of our system as 

This would correspond to the fact that we prepare our system in a equilibrium state at temperature 
^ for negative times and let the system evolve, for positive times, with the true Hamiltonian Ti 
which may be time dependent. If 'H(t) = Hi, then, the system will evolve in equilibrium and not 
otherwise. 



With eq. (^) in mind, we note that we can write 

U(T-i3,T) , , 

where T is assumed to be a large negative time (and not the temperature) and T — > — oo at the end. 
Using the semi-group properties of the time evolution operator, it is easy to see that the ensemble 
average of any operator can now be represented as 

{0)f, = Tr Pit) O 

Tr U{t, 0) U{T - if3, T) U{0, t) O 

Tr C/(r - T) 
T¥ U{T - if3, T) U{T, T') U{T\t) O U{t, T) 



Tr UiT - il3, T) U{T, T) U{T', T) 



(93) 



where we have introduced a large positive time T' and assume that T' — > oo at the end. This 
gives a nice representation to the ensemble average of any operator. Namely, we let the system 
evolve from a large negative time T to t where the appropriate operator O is inserted. The system 
then, evolves from t to a large positive time T' and back from T' to T and then, continues evolving 
along the imaginary branch from T to T — if3. Since the matrix elements of the time evolution 
operator can be given a path integral representation, it is clear that the ensemble average of any 
operator can also be given a path integral representation in this formalism corresponding to the 
specific contour in the complex time plane as described above. Although the specific contour has 
three branches - one along the real axis increasing with time, the second also along the real axis 
decreasing with time and the third along the negative imaginary axis - in the limit T ^ —oo and 
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T' ^ CO, it can be shown that the third branch gets decoupled from the other two (the factors in 
the propagators connecting such branches are asymptoticahy damped). Consequently, in this limit, 
we are effectively deahng with two branches leading to the name "closed time path formalism" [12]. 
In this contour, then, the time integration has to be thought of as 

dt= dt+- dt_ (94) 

Jc J ~oo J —CO 

where the relative negative sign arises because time is decreasing in the second branch of the time 
contour. 



The doubling of the degrees of freedom, in this formalism, is now clear. To have a path integral 
description, we must specify the fields on both the branches of the contour. Or, equivalently, we 
can use just the positive branch and double the field degrees of freedom. Namely, corresponding 
to every original field, say 0+, we must introduce a second field (j)- and remember that the action 
for the (/>_ fields must have a relative negative sign arising from eq. (|9j), namely, that time is 
decreasing along the second branch. 



Scalar Field Theory 



Just as an example, let us study next the self-interacting scalar field theory in some detail. The 



Lagrangian density is the same as in eq. (30), but following the earlier discussion, we should take 



the complete Lagrangian density for the system to be 

£ = £(0+)-£(0_) (95) 

where 

The Feynman propagator can again be determined for this theory and would have a 2 x 2 matrix 
structure because of the doubling of the field degrees of freedom. It can be determined subject to 
compatibility with the KMS conditions and has the form in the momentum space 

G{k) = ( ] (97) 



with 



G-+{k) G—{k) 



G++{k) = 5 -2i7mB{\k'^\)6{e -m^) 



G+-{k) = -2i7r(^e{-k^)+nB{\k°\))6{k^ -m 
G-+{k) = -2i7r(^eik^) + nB{\k°\))s{k'^ -m 



2^ 

2\ 



G—(k) = -i ^ 2mnBi\k^mk^ - m^) (98) 

^Z _ j^Z _ 
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There are several things to note from the structure of this propagator. First, as in the case 
of the propagator in thermo field dynamics, here, too, we see that the propagator naturally is a 
sum of two parts ~ the temperature independent part and the temperature dependent part. But, 
more interestingly, here the propagator has the simplification that the temperature dependent part 
of every component is the same which leads to various simplifications in actual studies of thermal 
quantities. Furthermore, not all the components of the propagator are independent. In fact, it is 
easily seen that (this can be traced back to their definition) 

G++ik) + G__(A;) = G+-{k) + G-+{k) 

These are known as the causal propagators of the theory and are useful in diagrammatic evaluation. 
There is, of course, another kind of propagator, conventionally known as the physical propagators 
and is defined as 

«(^-) - ( ) 

where Ga, Gr and Gc are known as the advanced, retarded and the correlated Greens functions. 
These are quite useful in the study of various phenomena such as the linear response theory. The 
important thing to observe is that the causal and the physical propagators are connected through 
a unitary transformation 

G{k)=QGQ^ (100) 

where 

It can be determined from this that, at the tree level, 

GA{k) ^ 
Gnik) 



J^2 _ ^2 _ ^g^O 
1 



^2 _ _|_ ^g^O 

Gc{k) = -2iTT(l + 2nB{\k^\))5{k'^ -m"^) (102) 



as they should be. 



The diagrammatic calculations can now be easily understood in this formalism. The vertices 
can be read out from the Lagrangian density in eq. (|9^). There are two kinds of vertices, one for 
the original fields, and the other for the doubled fields, </>_. The vertices for the 4>- fields are 
the same as those for the 4>+ fields except for a relative sign. With the vertices and the causal 
propagators, one can now carry out the calculation of any observable to any order in perturbation 
theory. As before, we note that, although there is no coupling between the and (p- fields at the 
tree level, higher order corrections would, in general, couple them. 

As an example, let us calculate the one loop mass correction in this theory. There will be two 
such diagrams to calculate - one for the (p+ field and the other for the 0_ field. The mass correction 
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for the 0+ field is readily seen to be 



\ f d^k / 1 



2mnB{\k^\)5{k^ - 



2 J (27r)4 VA;2 -m2 + ie 
where it is easily seen that the temperature independent part has the form 



i(Amg + Am|) (103) 



Amg = - / — 104 

4 7 2^ 3 cjfe 



while the temperature dependent part is given by 

^""^ " 2 7 (27r)3 a;^ "27 (27r)3 e/^-^ - 1 ^^"^^ 



These can be compared with the corresponding terms in eq. (35). We can also calculate the mass 
correction for the </>_ field. With a little bit of analysis, it is seen that 

Am^ = Am^ 



6 Feynman Parameterization 



So far, we have described the various formalisms that can be used to do calculations at finite 
temperature. However, actual calculations lead to many subtle, but interesting features of theories 
at finite temperature. One immediate and obvious feature, of course, is that finite temperature 
effects break Lorentz invariance. Namely, in studying a system at finite temperature, one has to 
go to a specific frame where the heat bath is at rest and, consequently, Lorentz invariance will no 
longer be manifest. This is, of course, already manifest at the level of propagators. For example, the 
structure of the propagators in eqs. (|8^) or ( p8| ) clearly displays a Lorentz non-invariant structure. 
The consequence of this is that an amplitude calculated at finite temperature, say for example, 
the self-energy Yi{p^ ^p) depends on the external energy and momentum independently. In fact, the 
self-energy becomes a non-analytic function of these two variables at the origin and two different 
ways of approaching the origin in this space leads to distinct plasmon and screening masses [23]. 
Thus, such non-analyticities are quite physical and their origin can be traced back to the fact 
that, at finite temperature, there are new channels of reactions possible leading to new branch 
cuts which give rise to such discontinuities [23, 10]. (To be absolutely fair, it is worth noting that 
statistical mechanics can be formulated in a covariant way. In such a case, one finds that there is a 
larger number of Lorentz invariant variables that can be defined on which amplitudes can depend. 
The non-analyticity in and p can then be translated to a non-analyticity in these new, Lorentz 
invariant variables [23].) 
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There are, of course, some other kinds of subtlety that arise which influence the calculations 
directly at finite temperature. We will discuss one such subtlety in this section. Let us note that 
a particularly useful formula in the evaluation of amplitudes at zero temperature is the Feynman 
combination formula given by 

1 1 _ /"i dx 

A + ieB + ie~ Jo [x{A + ie) + {1 - x){B + ie)]'^ ^ ' 

This can be directly checked by evaluating the x integral on the right hand side. 



This formula is extremely useful and works at zero temperature mainly because the Feynman 
propagators have the same analytic structure, namely, they have the same "ie" dependence. In 
contrast, we note that the finite temperature propagators contain delta functions (see eqs. (^) 
and (p8|)) and recalling that 

S(x) = lim -— 1 : ) 

e^0+ ZlTT \x — le x + le J 

we recognize that, at finite temperature, the propagators no longer have the same "ie" dependence. 
Consequently, in evaluating Feynman amplitudes at finite temperature, we have to combine de- 
nominators which do not necessarily have the same "ie" dependence. Keeping this in mind, let us 
examine the combination of two different denominators with arbitrary analytic dependence. 



Without loss of generality, let us choose a, (3 = ±1 and note that 

dx 1 1 



[x{A + iae) + {I - x){B + if3€)Y - B) + i{a - (3)e x{A + iae) + {I - x){B + ipej q 

~ A + iae B + i/Se ^^^^^ 

This is, of course, the standard Feynman combination formula. However, let us note that this will 
not hold if < xo < 1 such that 

xo = -^, (3A = aB (108) 
p — a 

because, in such a case, the integrand will have a singularity on the real j;-axis inside the interval 
of integration. In this case, we have 

dx dx 



[x{A + iae) + {l-x){B + ipe)f Jo [x{A - B + i{a - f3)e) + B + ifSef 



/3^-^ \ dx 



lim / + , I , 

ri->o\Jo J^+^J [x{A -B + i{a- (5)e) +B + if3ef 

A + iae B + il3e A-B + i{a-f3)e ^ ' 
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In other words, when the parameters of the integrand satisfy eq. ( |10^ ), the Feynman combination 
formula of eq. ( |106| ) will modify and the general formula follows from eq. (|109D to be [13] 

1 1 /•! dx (a-p)s(^f3A-aB) 



A + iae B + ifde Jo [x{A + iae) + (1 - x){B + if3e)Y A - B + i{a - P)e 



Note that, condition (|108| ) can only be satisfied (with < xq < 1) if a and P are of opposite sign. 
Indeed, let us note from eq. ( |11[1| ) that the second term vanishes when a = ^9 = 1 as is the case 
at zero temperature. Namely, when propagators with identical "ze" dependence are combined, the 



standard combination formula of eq. ( 106 ) holds. However, if denominators with opposite "ie" 
dependence are combined, the correct combination formula involves a second term. This is quite 
crucial at finite temperature and without this second term, one ends up with a wrong result as was 
discovered in finite temperature calculations the hard way [28]. 



7 Large Gauge Invariance 



Gauge theories are beautiful theories which describe physical forces in a natural manner and because 
of their rich structure, the study of gauge theories at finite temperature is quite interesting in itself. 
However, to avoid getting into technicalities, we will not discuss the intricacies of such theories 
either at zero temperature or at finite temperature. Rather, we will discuss a simple quantum 
mechanical model, in this section, to bring out some of the new features that temperature brings 
into such theories - features which are very different from what we expect at zero temperature. 

To motivate, let us note that gauge invariance is realized as an internal symmetry in quantum 
mechanical systems. Consequently, we do not expect a macroscopic external surrounding such 
as a heat bath to modify gauge invariance. This is more or less what is also found by explicit 
computations at finite temperature, namely, that gauge invariance and Ward identities continue 
to hold even at finite temperature [29]. This is certainly the case when one is talking about small 
gauge transformations for which the parameters of transformation vanish at infinity. 

However, there is a second class of gauge invariance, commonly known as large gauge invariance 
where the parameters do not vanish at infinity and this brings in some new topological character 
to physical theories. For example, let us consider a 2 + 1 dimensional Chern-Simons theory of the 
form 

£ = MCcS + -^^fermion 

= Me^"^ tr A^{d,Ax - ^A.A^) + ?(7^(i5^ - gA^) - m)i> (111) 

where M is a mass parameter, a matrix valued non-Abelian gauge field and "tr" stands for the 
matrix trace. The first term, on the right hand side, is known as the Chern-Simons term which 
exists only in odd space-time dimensions. We can, of course, also add a Maxwell like term to the 
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Lagrangian and, in that case, the Chern-Simons term behaves Uke a mass term for the gauge field. 
Consequently, such a term is also known as a topological mass term [30] (topological because it does 
not involve the metric). For simplicity of discussion, however, we will not include a Maxwell like 
term to the Lagrangian. 

Under a gauge transformation of the form 



^ U-^A^U-^U-^d^U (112) 



it is straightforward to check that the action in eq. (Ill) is not invariant, rather it changes as 



S= d^xC^S + ^^2mW (113) 
J 

where 

W = J d?x e^""^ tidf,UU-^duUU-^dxUU'^ (114) 

is known as the winding number. It is a topological quantity which is an integer (Basically, the 
fermion Lagrangian density is invariant under the gauge transformations, but the Chern-Simons 
term changes by a total divergence which does not vanish if the gauge transformations do not vanish 
at infinity. Consequently, the winding number counts the number of times the gauge transforma- 
tions wrap around the sphere.). For small gauge transformations, the winding number vanishes 
since the gauge transformations vanish at infinity. 



Let us note from eq. ( |113| ) that even though the action is not invariant under a large gauge 

2 

transformation, if M is quantized in units of the change in the action would be a multiple of 
2i7r and, consequently, the path integral would be invariant under a large gauge transformation. 
Thus, we have the constraint coming from the consistency of the theory that the coefficient of 
the Chern-Simons term must be quantized. We have derived this conclusion from an analysis of 
the tree level behavior of the theory and we have to worry if the quantum corrections can change 
the behavior of the theory. At zero temperature, an analysis of the quantum corrections shows 
that the theory continues to be well defined with the tree level quantization of the Chern-Simons 
coefficient provided the number of fermion flavors is even. The even number of fermion flavors is 
also necessary for a global anomaly of the theory to vanish and so, everything is well understood 
at zero temperature. 



At finite temperature, however, the situation appears to change drastically. Namely, the 
fermions induce a temperature dependent Chern-Simons term effectively making [31] 

mNf B\in\ 

M^M-- -4tanh-J— ^ 115) 

Att 2\m\ 2 ^ ^ 

Here, Nj is the number of fermion flavors and this shows that, at zero temperature (/3 — > oo), 
AI changes by an integer (in units of g'^/An) for an even number of flavors. However, at finite 
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temperature, this becomes a continuous function of temperature and, consequently, it is clear that 
it can no longer be an integer for arbitrary values of the temperature. It seem, therefore, that 
temperature would lead to a breaking of large gauge invariance in such a system. This is, on 
the other hand, completely counter intuitive considering that temperature should have no direct 
influence on gauge invariance of the theory. 



C-S Theory in + 1 Dimension 



As we have noted, Chern-Simons terms can exist in odd space-time dimensions. Consequently, let 
us try to understand this puzzle of large gauge invariance in a simple quantum mechanical theory. 
Let us consider a simple theory of an interacting massive fermion with a Chern-Simons term in 
0-1-1 dimension described by [14, 32] 

L = i^j{idt- A-m)i)j - kA (116) 

Here, j = 1, 2, • • • , Nf labels the fermion flavors. There are several things to note from this. First, 
we are considering an Abelian gauge field for simplicity. Second, in this simple model, the gauge 
field has no dynamics (in -|- 1 dimension the field strength is zero) and, therefore, we do not have 
to get into the intricacies of gauge theories. There is no Dirac matrix in -|- 1 dimension as well 
making the fermion part of the theory quite simple as well. And, finally, the Chern-Simons term, 
in this case, is a linear field so that we can, in fact, think of the gauge field as an auxiliary field. 



In spite of the simplicity of this theory, it displays a rich structure including all the properties 
of the 2 -|- 1 dimensional theory that we have discussed earlier. For example, let us note that under 
a gauge transformation 

^. ^ e-^^W^i, A^A + dtX{t) (117) 

the fermion part of the Lagrangian is invariant, but the Chern-Simons term changes by a total 
derivative giving 

S = [ dtL ^ S -2ttkN (118) 



where 

N 



l-JdtdtX{t) (119) 



is the winding number and is an integer which vanishes for small gauge transformations. Let us 
note that a large gauge transformation can have a parametric form of the form, say, 

X{t) = -iN log (^[^) (120) 

The fact that N has to be an integer can be easily seen to arise from the requirement of single- 
valuedness for the fermion field. Once again, in light of our earlier discussion, it is clear from eq. 
(|118|) that the theory is meaningful only if /t, the coefficient of the Chern-Simons term, is an integer. 
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Let us assume, for simplicity, that m > and compute the correction to the photon one-point 
function arising from the fermion loop at zero temperature. 

i,. = -H)iv, /f^ '<^-+;'', (121) 

This shows that, as a result of the quantum correction, the coefficient of the Chern-Simons term 
would change as 

Nf 

As in 2 + 1 dimensions, it is clear that the coefficient of the Chern-Simons term would continue to 
be quantized and large gauge invariance would hold if the number of fermion flavors is even. At 
zero temperature, we can also calculate the higher point functions due to the fermions in the theory 
and they all vanish. This has a simple explanation following from the small gauge invariance of the 
theory. Namely, suppose we had a nonzero two point function, then, it would imply a quadratic 
term in the effective action of the form 

r2 = ^ y dti dt2 A{ti)F{ti - t2)A{t2) (122) 

Furthermore, invariance under a small gauge transformation would imply 

ST2 = - J dtidt2X{ti)dt,F{ti-t2)A{t2) = (123) 

The solution to this equation is that F = so that there cannot be a quadratic term in the effective 
action which would be local and yet be invariant under small gauge transformations. A similar 
analysis would show that small gauge invariance does not allow any higher point function to exist 
at zero temperature. 



Let us also note that eq. ( |123|) has another solution, namely, 

-^'(^1 — ^2) = constant 

In such a case, however, the quadratic action becomes non-extensive, namely, it is the square of 
an action. We do not expect such terms to arise at zero temperature and hence the constant has 
to vanish for vanishing temperature. As we will see next, the constant does not have to vanish at 
finite temperature and we can have non-vanishing higher point functions implying a non-extensive 
structure of the effective action. 



The fermion propagator at finite temperature (in the real time formalism) has the form [10] 

S{p) = (p + m)(^ \ . 27mF{\p\)6{p'^ - in^)] 

VP — ?Ti^ + / 

1 

= 2TTnF(m)6(p - m) (124) 

p-m + ie ) yi' ) \ 1 

and the structure of the effective action can be studied in the momentum space in a straightforward 
manner. However, in this simple model, it is much easier to analyze the amplitudes in the coordinate 
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space. Let us note that the coordinate space structure of the fermion propagator is quite simple, 
namely, 

^^^^ " / ^ { p-m + ie ~ 2^^^M'5(P - m)J = {0{t) - nF{m))e-''^' (125) 
In fact, the calculation of the one point function is trivial now 

ill = -(-i)iV/5(0) = ^ tanh ^ (126) 

This shows that the behavior of this theory is completely parallel to the 2 + 1 dimensional theory 
in that, it would suggest 

K ^ K tanh 

2 2 

and it would appear that large gauge invariance would not hold at finite temperature. 

Let us next calculate the two point function at finite temperature. 

ih = -{-i?^S{ti-t2)S{t2-ti) 
= nFim){l — npifn)) 

= A,^,,^^ = ll.i2pi (127) 

8 2 2 2\ (3 dm ^ ^ 

This shows that the two point function is a constant as we had noted earlier implying that the 
quadratic term in the effective action would be non-extensive. 

Similarly, we can also calculate the three point function trivially and it has the form 

iNf , f3m , opm 11 d^{ili] 



ib = ^ tanh ^ sech^^ = Is (j) 

In fact, all the higher point functions can be worked out in a systematic manner. But, let us 
observe a simple method of computation for these. We note that because of the gauge invariance 
(Ward identity), the amplitudes cannot depend on the external time coordinates as is clear from 
the calculations of the lower point functions. Therefore, we can always simplify the calculation 
by choosing a particular time ordering convenient to us. Second, since we are evaluating a loop 
diagram (a fermion loop) the initial and the final time coordinates are the same and, consequently, 
the phase factors in the propagator ( [L25D drop out. Therefore, let us define a simplified propagator 
without the phase factor as 

S{t) = e{t) - npim) (129) 

so that we have 

S(t > 0) = 1 -nF(m), S{t <id) = -nF{m) (130) 
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Then, it is clear that with the choice of the time ordering, ti > t2, we can write 
dSih - 12) 



dm 
dS{t2 - h) 
dm 



(3S{ti - H)S{h - t2) ti>t2> h 
(3S{t2-H)S{H-ti) ti>t2>H (131) 



In other words, this shows that differentiation of a fermionic propagator with respect to the 
mass of the fermion is equivalent to introducing an external photon vertex (and, therefore, another 
fermion propagator as well) up to constants. This is the analogue of the Ward identity in QED 
in four dimensions except that it is much simpler. From this relation, it is clear that if we take a 
n-point function and differentiate this with respect to the fermion mass, then, that is equivalent 
to adding another external photon vertex in all possible positions. Namely, it should give us the 
(n + l)-point function up to constants. Working out the details, we have, 

^ = -^0{n + l)/.+i (132) 

Therefore, the (n + l)-point function is related to the n-point function recursively and, conse- 
quently, all the amplitudes are related to the one point function which we have already calculated. 
(Incidentally, this is already reflected in eqs. (|12^,[12^)). 



With this, we can now determine the full effective action of the theory at finite temperature to 

be 



-I 

n 



/ : — — tanh 

O /-^ rr,\ \ f)rrr, / 9 



where we have defined 



! \dm 

, a f3m a\ 

-iN f log ( cos 2 "I" ^ tanh sm - 1 (133) 



a = j dtA{t) (134) 



There are several things to note from this result. First of all, the higher point functions are 
no longer vanishing at finite temperature and give rise to a non-extensive structure of the effective 
action. More importantly, when we include all the higher point functions, the complete effective 
action is invariant under large gauge transformations, namely, under 

+ (135) 

the effective action changes as 

F ^ r + NNf-K (136) 
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which leaves the path integral invariant for an even number of fermion flavors. This clarifies the 
puzzle of large gauge invariance at finite temperature in this model. Namely, when we are talking 
about large changes (large gauge transformations), we cannot ignore higher order terms if they 
exist. This may provide a resolution to the large gauge invariance puzzle in the 2 + 1 dimensional 
theory as well. However, in spite of several nice analysis [33], this puzzle has not yet been settled 
in all its generality in the 2 + 1 dimensional case. 



Exact Result 



In the earlier section, we discussed a perturbative method of calculating the effective action at 
finite temperature which clarified the puzzle of large gauge invariance. However, this quantum 
mechanical model is simple enough that we can also evaluate the effective action directly and, 
therefore, it is worth asking how the perturbative calculations compare with the exact result. 

The exact evaluation of the effective action can be done easily using the imaginary time for- 
malism. But, first, let us note that the fermionic part of the Lagrangian in eq. ( [116D has the 
form 

Lf = lp{idt- A-m)^^ (137) 

where we have suppressed the fermion flavor index for simplicity. Let us note that if we make a 
field redefinition of the form 

ij{t) =e-'^o'^''^^''^it) (138) 
then, the fermionic part of the Lagrangian becomes free, namely, 

Lf = i>{idt-m)i) (139) 



This is a free theory and, therefore, the path integral can be easily evaluated. However, we have 
to remember that the field redefinition in ( |138 ) changes the periodicity condition for the fermion 



fields. Since the original fermion field was expected to satisfy anti-periodicity 

^{(3) = -V'(O) 
it follows now that the new fields must satisfy 

^(/3) = -e-*'^V'(0) (140) 



Consequently, the path integral for the free theory ( |139D has to be evaluated subject to the peri- 
odicity condition of ( 1401 ). 



Although the periodicity condition (140) appears to be complicated, it is well known that the 



effect can be absorbed by introducing a chemical potential [10], in the present case, of the form 
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With the addition of this chemical potential, the path integral can be evaluated subject to the 
usual anti-periodicity condition. The effective action can now be easily determined 



8 Supersymmetry Breaking 



One of the reasons for studying finite temperature field theory is to understand questions such as 
phase transitions in such systems. It is by now well understood that most field theoretic models 
of spontaneous symmetry breaking display a phase structure much like what one sees in a magnet, 
namely, above a certain critical temperature, the system is in a symmetric phase while below the 
critical temperature, the system is in a broken symmetry phase. Thus, temperature has the almost 
universal effect that if a symmetry is spontaneously broken at low temperature, it is restored 
at temperatures above a certain critical value. Qualitatively, it can be understood as follows. 
Temperature, particularly high temperature, provides a lot of thermal energy to a physical system 
to wash out any structure in the zero temperature potential which may be responsible for symmetry 
breaking. There is, however, one class of symmetries where temperature has the inverse effect, 
namely, in a super symmetric theory, a symmetric phase at low temperature goes to a broken phase 
at high temperature. (Of course, if supersymmetry is broken at low temperature, it continues to 
be broken even at high temperature.) We will discuss this phenomenon with a simple quantum 
mechanical model in this section. 



Supersymmetric Oscillator at T = 

Let us note that supersymmetry is an ultimate form of symmetry that one can dream of, namely, 
it transforms bosons into fermions and vice versa [34-35]. To introduce supersymmetry, let us 
consider a simple quantum mechanical model, commonly known as the supersymmetric oscillator 
[16]. It consists of a bosonic and a fermionic oscillator of the same frequency. Therefore, we can 
write the Hamiltonian, for the system as 



r 




(142) 



which coincides with the perturbative result of eq. (|134| ) . 




(143) 
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where ub and ap describe, respectively, the bosonic and the fermionic annihilation operators. 



The immediate thing to note from the structure of the Hamiltonian in eq. (|143| ) is that there is 
no zero point energy. We will see this shortly as a general feature of supersymmetric theories. Let 
us also define two fermionic operators of the form 



(144) 



With the usual canonical commutation relations for the bosonic operators (see eq. (|6q)) and anti- 
commutation relations for the fermionic operators (see eq. (|4^)), it is easy to check that 

[Q,H]=0 = \q,H 

Namely, these fermionic operators are conserved. In fact, together with the Hamiltonian, they 
satisfy the algebra (it is straightforward to check this) 



[Q,H] 

[Q,Q]+ 
Q,Q , 



: 
: 

u> 



Q,H 
Q,Q 



(145) 



Such an algebra, where both commutators and anti-commutators are involved (or alternately, where 
there is a grading of the multiplication rule of the algebra), is known as a graded Lie algebra and 
supersymmetric theories are realizations of graded Lie algebras. 



As we know from the study of symmetries, conserved quantities generate infinitesimal symme- 
tries of the theory. Since both Q and Q are conserved, it is worth asking what kind of symmetry 
transformations of the theory they generate. In fact, let us keep in mind that they are fermionic 
operators and hence the symmetry they will generate cannot be conventional. Explicitly, we can 
check that 



Q,aB 
[Q,aB 

Q,aB 
Q, ap 



= [Q,aF]_ 







Q,a\ 



a 



t 

F 
Qb 



(146) 



Namely, Q and Q take bosonic operators to fermionic ones and vice versa which is the bench 
mark of supersymmetry. Thus, our Hamiltonian in eq. ( |143D is invariant under supersymmetric 



transformations of the form (146). 
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There are several things to note from the structure of the supersymmetry algebra in eq. ( |145| ). 
First, the energy eigenvalues of our supersymmetric theory have to be positive semi-definite since 
the operator on the left hand side of the last relation in ( |145| ) is. Furthermore, if the ground state 
is supersymmetric satisfying 

Q|0) = = Q|0) (147) 

then, the ground state will have vanishing energy, as we had pointed out earlier as the case for our 
system. Both these results are, in fact, quite general for any supersymmetric theory. We also note 
from the structure of the algebra that the spectrum of the Hamiltonian will be doubly degenerate 
except for the ground state. Namely, if \tp) is an eigenstate of the Hamiltonian, then, Q\iIj) (or, 
Q\il)) - only one of them would be nontrivial depending on the form of \tp)) would be degenerate 
in energy. 



Let us, in fact, examine some of these general results explicitly. The spectrum of the Hamiltonian 
in eq. (|143| ) is, in fact, quite straightforward. The Hilbert space is a product space containing 
bosonic and fermionic oscillator states and a general state has the structure 



(4)"" (4) 



np 



|ns, np) = \nB) ® jn^) = ^ J-f ^ |0, 0) (148) 

with energy eigenvalues 

EnB,np = (^{riB + nF), = 0,1; ns = 0,1,2,- •• (149) 

where the ground state is expected to satisfy 

aB\0) = = aplO) (150) 
We note that an immediate consequence of ( |150| ) is that 

Q|0) = = Q|0) 

and, consequently, the ground state is supersymmetric and that the ground state energy is seen 
from (149) to vanish. All the higher states have positive energy. Furthermore, we note that all the 
states (except the ground state) of the form \nB,0) and \nB — 1,1) are degenerate in energy. Let 
us also note the effect of Q and Q acting on the states of the Hilbert space, namely. 



Q\nB,nF) 



^ynB + l\nB + l,nF - I) if np 
ifni;' = 

\nB — l,nF + l) if ub ^ Oor, np ^ 1 



Q\nB,np) = J - ^ ^ -f n 1 ^^^^^ 

I U if = Dor, n^? = 1 



Supersymmetric Oscillator at T 7^ 



Let us next analyze the supersymmetric oscillator at finite temperature in the formalism of thermo 
field dynamics. As we had noted earlier, this is the ideal setting to discuss questions such as symme- 
try breaking. Let us note, even before carrying out the calculations, that we expect supersymmetry 
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to be broken at finite temperature. Intuitively, this is quite clear. Namely, supersymmetry takes 
bosons to fermions and vice versa and, consequently, any boundary condition that distinguishes 
between the two would lead to a breaking of this symmetry. Temperature, in fact, introduces such 
a condition, namely, bosons and fermions behave differently at finite temperature (they obey dis- 
tinctly different statistics). However, what is not clear a priori is whether such a breaking would 
be explicit or spontaneous. 

To study the system at finite temperature within the framework of thermo field dynamics, let 
us look at the complete system, including the tilde oscillators, described by 

H = H — H = a;(a^aB + a^pap) — 0^(0^05 + cl^pClf) (152) 

The Hilbert space of the doubled system has the structure 

\nB,nF;fiB,nF) = \nB,nF) (8) \fiB,nF) (153) 

The thermal vacuum can now be defined (as discussed in section 3). Let us define 

G{6b, Of) = -i9B(.P){aBaB - o^b^b) ~ iOF{P){aFaF - ap^p) 
with (see eqs. (^)and (|72|) ) 

tan epiP) = ^ tanh 0b(^) (154) 

Then, the thermal vacuum can be defined as 

|0,/3) = e-^^(^S'^^) |0) (155) 
This also allows us to calculate the thermal operators in a straightforward manner. 

Let us note next that the expectation value of the Hamiltonian in the thermal vacuum is given 

by 

EoiP) = (0, (3\H\0, P) = (0, P\u{a)B<^B + a}pap)\Q, (5) 

= a;(sinh2 0s(/3) + sin2 0i.(/3)) = ^^— -32^ (156) 

This shows that the energy of the thermal vacuum is nonzero for any finite temperature signaling 
that supersymmetry is broken. Furthermore, let us note that 

Q\0,(3) = aiaF\0,P) = ^^^^^\nBiP) = l,nFiP)=0;nBiP) = 0,hpiP) = l) 

vl — e ^'^'^ 

_ , e-/3c^/2 

Q|0,/3) = a^paB\0,P) = ^===\nB{P) = 0,nF{P) = l;nB{P) = l,nF{P)=0) (157) 

vl — e ^^"^ 
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This, in fact, shows that supersymmetry breaking is spontaneous at finite temperature and the new 
states on the right hand side of (157) would correspond to the appropriate quasi particle Goldstino 
states associated with such a symmetry breaking. 



There are various other order parameters for the breaking of supersymmetry and all of them 
lead to the same conclusion that supersymmetry is spontaneously broken at finite temperature [16]. 



9 Conclusion 



In this article, we have tried to describe some of the interesting features of finite temperature field 
theories. There are, of course, many more topics that we have not been able to discuss. However, 
it is our hope that the topics discussed, in this article, would raise the curiosity of the readers to 
pursue various other questions in this field. 

This work was supported in part by the U.S. Dept. of Energy Grant DE-FG 02-91ER40685. 
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